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We report the results of searches for B + — > J/ipn' K + and B° — > J/ipn' Kg decays, using a sample 
of 388 xlO 6 BB pairs collected with the Belle detector at the T(4S) resonance. No statistically 
significant signal is found for either of the two decay modes and upper limits for the branching 
fractions are determined to be B{B + -» J/ipri'K + ) < 8.8 x 10" 5 and B(B° -> J/ipr)'K° s ) < 2.5 x 1(T 5 
at 90% confidence level. 

PACS numbers: 13.25.Hw,14.40.Gx,14.40.Nd 



Studies of exclusive B meson decays to charmonium 
play an important role in exploring CP violation [lj 
and in observations of new resonant states that include 
a (cc) pair The decay B — > J/vjrj'K proceeds 
via a Cabibbo-allowed and color-suppressed transition 
(b — > ccs) with (ss) quark popping. A branching frac- 
tion comparable to that of the decay B — > J/vj<f>K Q is 
therefore expected. 

The B —> J jvbrj ' K decay modes are of particular in- 
terest in the search for hybrid charmonium states, ip g , 
which are excited gluonic (ccg) states, which may decay 
to J J ifj + (?/'),</>) [5[. These hybrid states are expected 
to be produced in B decays such as B — » ip g K. A near- 
threshold u> J/vj enhancement in B — > KujJ/ip decays 
observed by Belle Q has some properties similar to those 
predicted for cc-gluon hybrid states 0. A state, F(4260), 
recently observed by BaBar in e + e~ — > F(4260)7(isr) 
transitions @, is also a candidate for such a ip g state 
where 7(isr) denotes an initial state radiation photon. 
An enhanced decay rate to J/iprj^ modes would be a 
supporting evidence for this assignment. 

In this paper, we report results on searches for the de- 
cay modes B+ -> J/ipri'K + and B Q -> J/iprj'K s [|| in a 
sample of 357 fb _1 containing 388 xlO 6 BB pairs accu- 
mulated at the T(4S I ) resonance with the Belle detector 
Q at the KEKB energy asymmetric e + e~ collider [if] ]. 
These are the first results ever presented for these decay 
modes. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector, a 50- 
layer central drift chamber (CDC), an array of aero- 
gel threshold Cerenkov counters (ACC), a barrel- like ar- 
rangement of timc-of- flight scintillation counters (TOF), 



and an electromagnetic calorimeter comprised of CsI(Tl) 
crystals (ECL). These detectors are located inside a su- 
perconducting solenoid coil that provides a 1.5 T mag- 
netic field. An iron flux-return located outside of the 
coil is instrumented to detect K® mesons and to identify 
muons. 

Candidates for these two decay modes are recon- 
structed with the decay chains J/ifi — > (/ = e,/i), 
T]' — > 7]7r + 7r _ , r\ — > 77 and Kg — > 7r + 7r~. 

Events are required to pass a basic hadronic event 
selection [11] . To suppress continuum background 
(e + e~ — > qq, where q — u,d,s,c), we require the ratio 
of the second to zeroth Fox- Wolfram moments 12] to be 
less than 0.5 and the absolute value of the cosine of the 
angle between the B meson and beam direction in the 
center-of-mass system to be less than 0.85. 

The selection criteria for the J/ip decaying to are 
identical to those used in our previous papers (llj . J/ijj 
candidates are pairs of oppositely charged tracks that 
originate from a region within 5 cm of the nominal in- 
teraction point (IP) along the beam direction and are 
positively identified as leptons. In order to reduce the 
effect of bremsstrahlung or final state radiation, photons 
detected in the ECL within 0.05 radians of the original 
e~ or e + direction are included in the calculation of the 
e + e~(-y) invariant mass. Because of the radiative low- 
mass tail, the J ftp candidates are required to be within an 
asymmetric invariant mass window: — 150(— 60) MeV/c 2 
< M e + e -( 7 )(M M+A( -)-m,j/ l/; <+36(+36) MeV/c 2 , where 
mj/^p is the nominal J /if) mass (l3j . In order to improve 
the momentum resolution of the J/tjj signal, a vertex and 
mass constrained fit to the reconstructed J/vj candidates 
is then performed and a loose cut on the vertex fit quality 
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is applied. 

In order to identify hadrons, a likelihood Li for each 
hadron type i (i — ir, K and p) is formed using informa- 
tion from the ACC, TOF and CDC (dE/dx). Charged 
tracks that were previously identified as electrons or 
muons are rejected in the hadron identification proce- 
dure. The kaon from the B meson and the pions from 77' 
are selected with the requirements of Lk/{Lk + L^) > 
0.4 and L 7r /(L 7r + Lk) > 0.1, which have efficiencies of 
82.5% and 89.4%, respectively. 

The 77 from 77' decay is reconstructed by combining 
two photons that do not match to any charged track. 
The calorimeter cluster energy, E 7 , is required to exceed 
50 MeV for both photons. To veto photons from it , 
we combine each photon candidate with another pho- 
ton (with E~ t > 50 MeV) and reject it if the invari- 
ant mass is consistent with a 7r°: |M 77 — m w o| < 18 
MeV/c 2 . For 77 candidates, cos0jj ol is used to suppress 
background, where #JJ c i is defined as the angle between 
the photon momentum in the 77 rest frame and the boost 
direction of 77 in the laboratory frame. Signal events have 
a uniform cos 0? , distribution, whereas random photon 
pairs have a distribution that peaks near ±1. We require 
IcosflJJJ < 0.94 (0.97) for charged (neutral) B decay. 
The invariant mass of the 77 candidate is required to be 
within 0.496 MeV/c 2 < M 77 < 0.582 MeV/c 2 (3cr of the 
nominal 77 mass). To improve the momentum resolution 
of the 77, we apply a mass constrained fit to the recon- 
structed 77. 

We reconstruct 77' candidates by combining the 77 and 
selected tt + it~ candidates. We require the invariant mass 
to be 0.940 MeV/c 2 < AV +7r - < 0.975 MeV/c 2 . To 
improve the momentum resolution, we apply a vertex 
and mass constrained fit on the 77' [14j |. 

For Kg candidates, we impose momentum-dependent 
requirements on the impact parameters from the IP 
for both Kg daughter tracks, the distance between the 
daughter tracks along the beam axis at the K s vertex, 
the difference of azimuthal angles between the K s mo- 
mentum and the direction of the K s vertex from the IP, 
and the flight length of the Kg. The invariant mass of the 
K s candidate is required to be within 16 MeV/c 2 (3c) of 
the K s mass. We apply vertex and mass constrained fits 
for the K s candidates to improve the momentum resolu- 
tion. 

These criteria maximize iV s i g / yj iV s i g + iVbkg , where 
iVsig is the number of expected signal events from sig- 
nal Monte Carlo (MC) samples with assumed branching 
fractions of 1.0 x 10~ 5 for both B + — > J/ipr]'K + and 
B° — > J/ipr)'K s modes, and A^kg is the number of ex- 
pected background events estimated from a MC sample 
of i?_B-pairs with a J/tp in the final state. 

B mesons are reconstructed by combining a J/tp, an 
77' and a K candidate. We identify B candidates us- 
ing two widely used kinematic variables calculated in 
the center-of-mass system: the beam-energy constrained 
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beam 



Pg) and the energy difference 
is the beam en- 



mass (Mbc 

(AE ee E B - Bbeam), where E hc 
ergy, Pg and Eb are the reconstructed momentum and 
energy of the B candidate. We select both B + and 
B° candidates within the range \AE\ < 0.20 GeV and 
Mbc > 5.20 GeV/c 2 for the final analysis. The signal 
regions are defined to be 5.27 GeV/c 2 < M hc < 5.29 
GeV/c 2 and \AE\ < 0.03 GeV for both B+ -> J/ipr]'K + 
and B° — > J/ipr/'Kg modes, which corresponds to three 
standard deviations based on the MC simulation. 

After all selection requirements, about 26.5% of the 
B + — > J/i/ji]'K + candidates have more than one entry 
per event; this occurs for 30.0% of the B° — > J/tprj'Kg 
candidates. For these events, the B candidate with the 
best vertex fit quality is used. From a study of MC sim- 
ulated events, we conclude that around 81.5% of multi- 
ple B + — > J/i(jt]'K + candidates and 70.0% of multiple 
B° — > J/ipri'Kg candidates are selected correctly by this 
procedure. 

The signal yields are extracted by maximizing the two 
dimensional extended likelihood function, 



V! 



n 



where N is the total number of candidate events, i is the 
identifier of the i-th event, Nk and Pk are the yield and 
probability density function (PDF) of the component fc, 
which corresponds to the signal and background. 

The signal PDFs for the two decay modes are mod- 
eled using a sum of two Gaussians for M^ c and a sum of 
two Gaussians plus a bifurcated Gaussian that describes 
the tail of the distribution in AE. The PDF parame- 
ters are initially determined using signal MC and subse- 
quently the primary Gaussian parameters are corrected 
using a control data sample of B + — » J/ipK* + decays 
with K* + — ► ir°K + . The signal shape parameters are 
kept fixed in the fits to the data. 

The dominant background comes from random combi- 
nation of J/tp, 77', and K candidates in BB events. The 
background from continuum is found to be small (a few 
%). A threshold function is used for the M hc PDF. 
For AE, we use a first-order polynomial function. The 
MC study shows that the background from the decays 
with similar topology that would make a peak in Mb c or 
AE distributions is negligible. 

In the fit, the value of Nk and the parameters for 
the combinatoric background are allowed to float. Fig- 
ures [l] and [2] show the (Mb C) AE) scatterplots and their 
projections for candidates after all selections are ap- 
plied. The fit results are superimposed on the projec- 
tions. There are 37 candidate events in the signal region 
for B+ -> J/ij)r]'K + and two for B° -> J/ipr/K%. 

Table [J summarizes the maximum-likelihood fit results 
for the signal (Y) and signal-region background (b) yields 
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FIG. 1: (a) The (Mbc, AE) scatterplot of B + -> J/ipr]'K + candidates and its projections onto (b) AE with 5.27 GeV/c 2 
< M bc < 5.29 GeV/c 2 and (c) M bc with |AE| < 0.03 GeV. The dashed lines and solid boxes indicate the signal regions. The 
curves are the result of the fit as described in the text. 



> 02 
a 

g.0.15 
HI 

< 0.1 



0.05 


-0.05 
-0.1 
-0.15 

•°i2 



: . J- .• N 



T....*» 



..*I.*....T..*J! 



••• 



• mm ■ * •= • 



5.22 5.24 5.26 5.28 5.3 
M bc (GeV/c 2 ) 




1.2-0.15-0.1-0.05 



0.05 0.1 0.15 0.2 
AE(GeV) 



- all 

signal j 
background 



(C> 



5.22 5.24 



5.26 5.28 5.3 
M bc (GeV/c*) 



FIG. 2: (a) The (M bc , AE) scatterplot of B° -> Jj^r(K% candidates and its projections onto (b) AE with 5.27 GeV/c 2 
< M bc < 5.29 GeV/c 2 and (c) M bc with \AE\ < 0.03 GeV. The dashed lines and the solid box indicate the signal regions. The 
curves are the result of the fit. 



TABLE I: Summary of the results. Y and b are the signal and expected total background yields in the signal region, Sig. is the 
statistical significance, no is the observed number of candidate events in the signal region, Y90 is the upper limit on the signal 
yield at 90% confidence level, e (error includes systematic error) is the detection efficiency and B is the branching fraction. 





Mode 




b 


Sig. 


n 


Y 90 


<%) 


B 


B+ 
B° 


-> J/ipr,'K + 


6-0±Slg 


22.0 ±0.7 
4.2 ±0.3 


1.8cr 

0.9cr 


37 
2 


28.1 
3.7 


3.9 ±0.6 

9 y+0.6 
z - ' -0.7 


< 8.8 x 10" b 

< 2.5 x 10" 5 



and their statistical errors. The statistical significance is 
defined as yf— 21n(£o/£ m ax), where £ m ax and £0 denote 
the maximum likelihood with the fitted signal yield and 
with the yield fixed at zero, respectively. 

No significant signal is found for either the B + — > 
J/4>'H I K + or B° — > J/iprj'Kg decay mode. We obtain 
upper limits on the yield at 90% confidence level (Y90) 
from number of observed candidates in the signal region 
(no) and b using the Feldman-Cousins method [lij . We 
account for systematic uncertainty due to uncertainties 
in the signal reconstruction efficiency and background es- 
timate using the method of Ref. [17| • 



The branching fraction is determined with Ns/[t x 
N BS x B(J/tp l+l~) x B(i -> r)ir+ir-) x B(r) -> 77)] 
for B+ -+ J/ipr]'K + and N s /[e x N BB x B(J/ip -» 

l+l~) X S(?7' nTT+TT-) X B(n 77) X B(2Tg 7T+7T-)] 

for B° J/ipr]'K%. Here N s is the signal yield, N BB 
is the number of BB pairs. We use the world aver- 
ages for the branching fractions of B(J/tp — * l + l~), 
B{q' -> nw+TT-), B(r) -> 77) and S(i?fg -> tt+tt-). The 
efficiencies (e) are determined from the signal MC sample 
with the same selection as used in the data. A three-body 
phase space model is employed for all three decay modes. 
The fractions of neutral and charged B mesons produced 
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TABLE II: Summary of systematic uncertainties (%) of the denominator in branching fraction calculation. 



Source 


B + -> J/*i/)r)'K + 


B" - J/HK S 


Tracking error 


±7.4 


±8.9 


PID (tt and K) 


±2.7 


±1.4 


Lepton ID 


±3.8 


±3.8 


7T° veto 


±3.0 


±3.0 


77 reconstruction 


±4.3 


±4.3 


Kg reconstruction 




±3.7 


Branching fractions 


±3.3 


±3.3 


MC statistics 


±0.2 


±0.3 


3-body decay model 


+9.4/-11.3 


+17.5/-22.8 


N BS 


±1.4 


±1.4 


Total 


+14.3/-15.6 


+21.4/-25.9 



in T(45) decays are assumed to be equal. 

The sources and sizes of systematic uncertainties in 
branching fraction calculation are summarized in Ta- 
ble [III The dominant sources are uncertainties in the 
three-body decay model, tracking efficiency and particle 
identification. For the error due to uncertainty in de- 
cay modeling for B+ -> J/tJjt]'K+ (B° -> J/~iprfK s ), the 
distribution in phase space is unknown. We conserva- 
tively assign the maximum variation of efficiency among 
the slices of M(J/ip,r)'), M(J/<ip,K-) and M(ry',i^") 



rea); KBN (contract No. 2P03B 01324, Poland); MIST 
(Russia); MHEST (Slovenia); SNSF (Switzerland); NSC 
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[M(J/V>,f/), M (J/i>, K s) and M(r/,K%)] as a system- 
atic uncertainty. The uncertainties in the tracking effi- 
ciency are estimated by linearly adding the momentum- 
dependent single track systematic errors. We assign un- 
certainties of about 1.3% per kaon, about 0.7% per pion, 
and about 1.9% per lepton for the particle identification. 
The r\ reconstruction uncertainty is determined by mea- 
suring the efficiency ratio between data and MC sam- 
ple in the inclusive 77 sample. The Kg reconstruction 
is checked by comparing the ratio of D + — > K s tt + and 
D + — > K~iy + tt~ yields. 

The systematic errors for the background yield are 
evaluated by varying each of the PDF parameters by its 
statistical error from the fit. 

In summary, we searched for B + — > J/ipi]'K + and 
B° — > J/ipj]'K s decays. No statistically significant signal 
was found for cither of the two decay modes; upper limits 
for the branching fractions are determined to be B(B + — > 
J/ipri'K + ) < 8.8 x 10~ 5 and B(B° -> J/^rj'K s ) < 2.5 x 
10~ 5 at 90% confidence level. 
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